Exotic plasmon modes of double layer graphene 
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We derive the plasmon dispersion in doped double layer graphene (DLG), made of two parallel 
graphene mono-layers with carrier densities m and 712 respectively, and an inter-layer separation 
of d. The linear chiral gapless single particle energy dispersion of graphene leads to exotic DLG 
plasmon properties with several unexpected experimentally observable characteristic features such 
as a nontrivial influence of an undoped (712 = 0) layer on the DLG plasmon dispersion and a strange 
influence of the second layer even in the weak-coupling d — > limit. Our predicted exotic DLG 
plasmon properties clearly distinguish graphene from the extensively studied usual parabolic 2D 
electron systems. 
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Graphene, a two dimensional (2D) layer of carbon 
atoms arranged in a honeycomb lattice, has attracted 
a great deal of attention because of its unique electronic 
properties arising from its chiral, gapless, linear band dis- 
persion Of particular importance and fundamental 
interest are those graphene properties which are qualita- 
tively different from the behavior of electrons (or holes) 
in the extensively studied 2D semiconductor systems. 
Some of these unique graphene properties [l| are its 'half- 
integral' integer quantized Hall effect, a carrier mobility 
independent of carrier density away from the charge neu- 
trality point (i.e. the Dirac point), and the finite min- 
imum conductivity around the Dirac point without any 
obvious metal-insulator transition. In this Letter, we pre- 
dict a new exotic collective property of graphene elec- 
trons, which distinguishes graphene qualitatively from 
the semiconductor-based parabolic 2D electron systems. 
We find the surprising result that the plasma modes of 
an interacting double layer graphene (DLG) system are 
completely different from the well studied double layer 
semiconductor quantum well plasmons [2]. Our finding 
is unexpected and interesting because the collective plas- 
mon excitations arise from the long-range Coulomb in- 
teraction which is the same in graphene and parabolic 
2D systems. The exotic DLG palsmon properties are ex- 
perimentally observable and should have important con- 
sequences for the many-body effects in doped multi-layer 
graphene structures. 

Plasma modes of 2D double layer structures have been 
extensively studied ever since the existence of an un- 
damped acoustic plasmon mode was predicted Q in semi- 
conductor double quantum well systems. It is known 
that when the two 2D layers are put in close proxim- 
ity with a high insulating barrier between them to pre- 
vent interlayer tunneling, the 2D plasmons are coupled by 
the interlayer Coulomb interaction leading to the forma- 
tion [2] of two branches of longitudinal collective excita- 
tion spectra called the optical plasmon (OP, cj + ~ g 1 / 2 , 



mon (AP, co- ~ <?), where the density fluctuations in 
each component oscillate in phase (OP) and out of phase 
(AP), respectively, relative to each other, assuming the 
two charge components to be the same. These collective 
modes of the double layer structures, which have been di- 
rectly experimentally observed [3|, play important roles 
in the many body properties such as screening and drag 

an. 

Here we consider the experimentally relevant issue of 
the collective mode dispersion in DLG formed by two par- 
allel single layer graphene (SLG) separated by a distance 
d. The DLG is fundamentally different from well studied 
bilayer graphene p because there is no inter-layer tunnel- 
ing, only inter-layer Coulomb interaction. Spatially sep- 
arated two-component DLG can be fabricated by folding 
an SLG over a high insulating substrate [7[. 

We calculate the DLG plasmon modes and their loss 
functions (spectral strength). We find many intriguing 
and unexpected features of coupled plasmon modes and 
their Landau damping in DLG. Even though we recover 
in the long wavelength (q — » 0) limit the well-known opti- 
cal (cj + ~ q 1 / 2 ) and acoustic (uj- ~ q) plasmons, the den- 
sity dependence is given by w\(q — > 0) ex ^/rTf + y^nj and 
(q — > 0) oc ^Jn^wij (y/n\ + ^/nj) compared to uo\ ex N 
and uPl (X n\ri2/N in the ordinary 2D system, where 
N = m + n2. When the interlayer Coulomb coupling is 
weak (i.e. for large separations, kpd ^> 1) the undamped 
DLG modes (oo±) become degenerate and have the same 
frequency as the uncoupled SLG plasmon [8], but inside 
the interband Landau damping regime this degenerate 
mode splits into two coupled modes. In the strong cou- 
pling regime (kpd <C 1) the cj_ mode approaches the 
intraband electron-hole continuum (i.e. ou- — > v^q) and 
is overdamped, but the o; + mode is shifted to higher fre- 
quency and in the long wavelength limit becomes the 
SLG plasmon with a density N = n\ + ri2 + 2^/77477,2 
instead of the density N = m + 772 as happens in reg- 
ular 2D double layer systems. We also find that the 



where q is the 2D wave vector) and the acoustic plas- DLG plasmon modes are heavily damped even in the 
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long wave length limit when the two layers have unequal 
density, even when one layer is undoped (i.e. n 2 = 0). 
We believe that our predictions may be easily observable 
via inelastic light scattering spectroscopy frequency- 
domain far infrared (or microwave) spectroscopy [lOj . or 
inelastic electron scattering spectroscopy [111, Elj . 

The plasmon modes are given by the poles of the 
density-density correlation function, or equivalently by 
the zeros of the dynamical dielectric function. For a dou- 
ble layer system the collective modes are given by the ze- 
ros of the generalized dielectric tensor ei m where I, m=l 
or 2 with 1,2 denoting the two layers. [2] The dielectric 
function ei m (q,uj) = Si m — ^ m II m is obtained within the 
mean field random phase approximation (RPA) in our 
theory, and the two component determinantal equation 
becomes 

e(q,uj) = [1 -^n(g)IIi(g,a;)][l -^2 2 (g)n 2 (g,a;)] 

- v 12 (q)v2i(q)^-i(q^)U 2 (q,u;). (1) 

where vu(q) and vi m (q) are respectively the intralayer 
and inter layer Coulomb interaction matrix elements. In 
our model, vu(q) = 27re 2 /(^); vi m (q) = vu(q)exp(-qd), 
with n as the background lattice dielectric constant. Fi- 
nally, Hi(q,uo) in Eq. (pQ) is the noninteracting SLG po- 
larizability function given by [8]: 



fsk — fs'k' 



L 2 1 — ' U) ■ 

kss 7 



IT] 



-F M ,(k,k'), (2) 



where g = g s g v (with g s = 2 and g v = 2 being spin 
and valley degeneracy), k' = k + q, s,s' = ±1 denote 
the band indices, e s k = svp|k| (vp being the Fermi ve- 
locity of graphene and h = 1 throughout this paper), 
and F ss /(k, k') is the overlap of states and given by 
F ss /(k, k') = |(1 + ss f cos 6), where 6 is the angle be- 
tween k and k', and f s k is the Fermi distribution func- 
tion, f sk = [exp{/3(e sk - fn)} + l]-\ with = l/k B T 
and \i\ the chemical potential of /-the layer. In the long 
wavelength limit (q — > 0) the polarizability becomes [8[ 
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where {Ef 1 = vf^Fi) is the Fermi wave vector (Fermi 
energy) of the l-th layer. 

The analytical formula for the long-wavelength plas- 
mon dispersion can be obtained in the general situation 
where the two layer have different Fermi wave vectors 
and Fermi energies by virtue of having different densities 
774, ^2- From Eqs. (pQ) and ([3|) it is possible to obtain the 
following long wavelength plasma modes of the coupled 
DLG: 
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FIG. 1: Calculated plasmon mode dispersions of coupled DLG 
for a fixed density and several layer separations. The upper 
(lower) solid line indicates the in-phase mode o;+ (the out of 
phase mode co-). Here we use the parameters: n\ — — 
10 12 cm- 2 and (a) d = 20 A (k F d = 0.35), (b) d = 100 
A (k F d = 1.8), (c) d = 300 A (k F d = 5.3), (d) d = 500 A 
(kpd = 8.9). The dot-dashed line indicates the plasmon mode 
dispersion of SLG with the same density. The dashed lines 
indicate the boundaries of SPEs (Landau damping region for 
intra- and interband electron- hole excitations). 



where r s = e 2 j kvf is the effective graphene fine structure 
constant. Thus, one recovers in a straightforward fashion 
the well-known optical (oj+ ~ q 1 ^ 2 ) and acoustic (co- ~ 
q) plasmons of a double layer system. In addition we have 
the density dependence of DLG plasmons as oo\(q — > 
0) oc y/ri\ + ^fni and uo 2 _(q — > 0) oc v /nin 2 /(v / nT + v^)' 
which are different from the density dependence of the 
corresponding plasmon dispersion in ordinary 2D systems 
[2[ where the uo 2 ^ depends on the total 2D electron density 
N = m + n2, (i.e. uo\ ex TV), and uj 2 _ oc ninz/N. 

In Fig. [T] we show the calculated DLG plasmon disper- 
sions (oo±) for equal densities of n\ = n 2 = 10 12 cm~ 2 
and several layer separations. The plasmon mode dis- 
persion of the uncoupled SLG with the same density is 
shown as the dot-dashed line. In Fig.[TJ we also show the 
electron-hole continua or single-particle excitation (SPE) 
region, which determines the decay (Landau damping) of 
the plasmon at a given frequency and wave vector. If the 
collective mode lies inside the SPE continua, we expect 
the mode to be Landau damped. The SPE continuum 
is defined by the nonzero value of the imaginary part 
of the total dielectric function, Im[e(g, a;)] ^ 0, because 
the damping of the plasmon by emitting an electron-hole 
single pair excitation is allowed in the region of nonzero 
lm[II(g,cj)]. For graphene, both intraband and inter- 
band SPE transitions are possible [8[ , and the boundaries 
(dashed lines) are given in Fig. [lj Due to the phase- 
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space restriction, the interband SPE continuum opens 
a gap in the long wavelength region satisfying the fol- 
lowing condition: vpq < uo < 2Ep — Vpq, in which the 
undamped plasmon modes exist. Since the normalized 
plasmon mode dispersion in Ep scales with normalized 
Fermi wave vector, our calculated results of Fig. [T] apply 
also for all other densities. 

As shown in Fig. ffl(a), f° r a sma H spatial separa- 
tion (kpd < 1), the frequency of the acoustic mode 
co- decreases compared with the uncoupled mode and 
approaches the upper boundary of intraband SPE (i.e. 
uo = vfq) while the optical mode uo+ shifts to higher 
energy. Especially, as d — > 0, the uo- mode becomes 
degenerate with the electron-hole continuum and loses 
its identity (i.e. uo- — > vpq). On the other hand, 
the uo+ mode becomes the SLG plasmon with a density 
N = rii + ^2 + 2^JniTi2 in the long wavelength limit. 
This is very interesting because, when the two layers with 
densities ni,ni become one combined system (d — > 0), 
the optical plasmon mode (a;+) should correspond to the 
SLG plasmon of total density N = n\ + ri2, which is pre- 
cisely what happens in regular 2D double layer system. 
As the spatial separation increases the optical (acoustic) 
plasmon frequency decreases (increases) slowly and the 
mode coupling occurs at the boundary of interband SPE 
(see Fig. OJb)). As d increases further the two modes 
become degenerate below the SPEi nte r (i.e. in the un- 
damped region), but the degenerate mode is divided into 
two coupled modes inside interband SPE (see Fig. []Jc) 
and (d)). For infinite separation, the two modes become 
uncoupled, and have the dispersion of SLG plasmon with 
the density n = ni = 712. 

In Fig. [2] we show the calculated DLG loss function (i.e. 
— Im[e(g, uo) -1 )) plotted in q-uo space. The loss function 
is related to the dynamical structure factor S(q,uo) by 

S(q,w) oc — Im e(q^uo)~ 1 . The dynamical structure fac- 
tor gives a direct measure of the spectral strength of the 
various elementary excitations. Thus, our calculated loss 
function can be measured in experiments such as inelas- 
tic electron and Raman scattering spectroscopies. When 
both Re[e] and Im[e] become zero (i.e. e(q,u) = 0, which 
defines the plasmon mode), the imaginary part of the 
inverse dielectric function, Im[e(g, ^) _1 ], is a ^-function 
with the strength W(q) = 7r[|<9Re[e(g, uo)]/ doo^^q)]' 1 , 
where uo p (q) is the plasmon dispersion. Thus the un- 
damped plasmon shows up as well defined ^-function 
peak in the loss function. The damped plasmon, how- 
ever, corresponds to a broadened peak in the loss func- 
tion - for larger broadening, the plasmon is overdamped, 
and there is no peak in the loss function. 

In Fig. [2] we show the — Im[e(<2, oo)~ x ] in density plots 
where darkness represents the mode spectral strength. 
The solid lines indicate the ^-function peaks and cor- 
respond to the well-defined undamped plasmon modes. 
The intraband electron-hole SPE continuum shows up 
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FIG. 2: The density plots of DLG loss function 
(— Im[e(<7, uo) -1 ]) in (q,uo) space for fixed densities of m — 
ri2 = 10 12 cra -2 and several layer separations d — 20, 100, 
300, 500A. The dashed lines indicate the boundaries of SPEs 
and the solid lines represent the undamped plasmon modes 
(^-function peaks in the loss function). 



as weak broad (incoherent) structure in Fig. [2] and car- 
ries small spectral weight. For a small spatial separa- 
tion only the in-phase mode oo+ carries any significant 
spectral strength. As d increases the out-of-phase mode 
carries substantial spectral weight even at long wave- 
lengths (q — > 0). Both oo± modes in general carry fi- 
nite spectral wei ghts and should be observable in exper- 
iments 

sua Etna. 

Since both co±(q) are greater than 
uo = VFq for all q, which is the upper boundary of the 
intraband electron-hole pair continuum, we expect these 
modes not to decay to intraband electron- hole pairs. But 
the modes enter into interband SPE at finite g c , which is 
given by the condition, oo±(q c ) = 2Ep — Vpq c . If q > q c , 
c(q,oo±) has a finite imaginary part and the plasmon 
modes become damped consequently (Landau damping). 
The plasmon mode inside the Landau damping region 
decays by emitting interband electron-hole pairs which 
is now allowed by energy- momentum conservation. Near 
the boundary of the interband SPE the mode damping 
is not significant. 

In Fig. [3] we show the calculated DLG plasmon 
mode dispersion and the corresponding loss function 
(— Im[e(<2, cj) -1 ]) for unequal layer densities at fixed layer 
separation (d = 100 A). As the layer densities are differ- 
ent the undamped plasmon region becomes smaller in q-uo 
space. As q increases the modes first enter SPE? nter at 
smaller q c and decay by producing interband electron- 
hole pairs in the layer 2. As g increases further the 
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FIG. 3: Calculated plasmon mode dispersions and 
— Im[e(g, cj) _1 ] of DLG for fixed layer separation (d — 100A) 
and for different layer densities. Here we use m — 10 12 cm~ 2 
and in (a)(c) 712/ n\ = 0.2 and in (b)(d) 712/ni = 0.5. 



FIG. 4: (a) Calculated plasmon mode dispersions of DLG for 
different layer separations for 712/ f ni =0 (i.e. 712 = and m = 
10 12 cra -2 ), and corresponding loss functions (— Im[e(g, u;) -1 ]) 
for (b) d = OA, (c) d = 10A, and (d) d = IOOA. 



modes enter into SPE^ nter and decay by exciting inter- 
band electron-hole pairs in both layers. As the density 
imbalance increases (i.e. n^l 'rt\ decreases) the acoustic 
mode approaches the boundary of intarband SPE 

and loses its spectral strength severely. But the in-phase 
mode 00+ remains a well defined peak inside SPE? nter . 

In Fig. ffl(a) we show the DLG plasmon modes in the 
case of ri2/ni = 0, i.e., the layer 2 is undoped (712 = 0) 
and the layer 1 has a finite density rt\ = 10 12 era -2 . 
Even though there are no free carriers in layer 2 we have 
two coupled plasmon modes (i.e we find two zeros in 
Re[e(g, oj)])\ one is above SPEi ntra and the other 

(00-) is degenerate with the boundary of SPEi ntra (i.e. 
uj(q) — vpq). More interestingly, when the layer sepa- 
ration is finite (i.e. d 7^ 0) the two plasmons merge at 
a finite wave vector q m (> kp). We find that the mode 
dispersion and q m are unchanged when d > 100 A. As 
d — > the dispersion of oo+ becomes exactly that of SLG 
plasmon with the same density, and the two modes oo± 
never meet. In Fig. Sfb)-(d) we show the loss function 
of DLG corresponding to the Fig. Ufa). Since all modes 
are inside of SPE 2 nter there is no undamped Landau re- 
gion, so the modes are damped by producing interband 
electron-hole pairs in layer 2 even in the long wavelength 
limit. The loss function shows that 00 - does not carry 
any spectral weight, and all spectral weight is carried by 
00+. For d = 0, the 00+ mode becomes a well defined 
peak only in the long wavelength limit and appears as a 
very broad peak for high wave vectors. As the separa- 



tion increases, 00+ is well defined below q m even inside 
SPEi nter and manifests well defined spectral peaks. Note 
that there is no significant spectral weight for q > q m . 

In conclusion, we investigate theoretically the dynam- 
ical response of DLG. We find interesting novel plasmon 
features unique to DLG which do not exist in ordinary 
2D double quantum well structures. 
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